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Abstract 



We study the sharp order of integrability of the exit position of Brow- 
nian motion from the planar domains V a = {(x,y) £ M x 1: a; > 
0, \y\ < Ax a }, < a < 1. Together with some simple good-A type 
arguments, this implies the order of integrability for the exit time of 
these domains; a result first proved for a — 1/2 by Banuelos, DeBlassie 
and Smits pQ and for general a by Li |5] . A sharp version of this result 
is also proved in higher dimensions. 
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1 Introduction 



For 7 with < 7 < tt, we denote by T 7 the right circular cone in R n of angle 
7. We let {Bt : t > 0} be n-dimensional Brownian motion and denote by 
E x and P x the expectation and the probability associated with this motion 
starting at x. We write T 7 = inf{t > : Bt T 7 }, so that T 7 is the first exit 
time of the Brownian motion from T 7 . In 1977, Burkholder [3] found the 
sharp order of integrability of T 7 . More precisely, he showed the existence 
of a critical constant p(pf, n) (given explicitly in terms of zeros of certain 
hypergeometric functions) such that, for z £ Tj, 

(1.1) E Z T^ 2 < 00, 

if and only if p < p{pi, n). In particular, for a cone in two dimensions (a case 
that had already been solved in 11) E Z T% /2 < 00 if and only if p < tt/ (27). 
Notice that by making the angle of the cone arbitrarily small we can make 
p arbitrarily large. 

The cone in two dimensions can be thought of as the domain above 
the graph of the function y = a\x\. It is then natural to study the order 
of integrability of exit times from other unbounded regions, in particular 
parabolas. Since any parabola is contained in a cone of arbitrarily small 
angle, the exit time for a parabola has finite moments of all orders. On the 
other hand, by comparing with rectangles, it is also easy to show that the 
exit time is not exponentially integrable. Bahuelos, DeBlassie and Smits 
showed in pQ that if T-p is the exit time of the Brownian motion from the 
parabola V = {(x,y): x > 0, \y\ < A\fx\ and z £ V, then there exist 
positive constants A\ and A<i such that 

(1.2) -Ax < liminft~i log [P z {r v > t}} 

t— too 

< lim sup log [P z {t-p > t}} < -A 2 . 

Thus if a < 1/3, then £'^[exp(6rp)] < 00 for each b > and, if a > 1/3, 
then E z [exp(bTp)] = 00 for each b > 0. This result was extended to all 
dimensions and to other unbounded regions by W. Li [Sj. More recently, 
M. Lifshits and Z. Shi QUI found the limit -I of log [P z {t v > t}], as 

1 /3 

t — > 00. Then, if b < I then ^[exp(6r^ )] < oo and if b > I then 

1/3 1/3 
E z {exp(bT-p )] = 00, (it is not known whether or not E z [exp(lTj, )] is finite 

or infinite) , so that Lifshits and Shi determine the sharp order of integrabil- 
ity of the exit time T-p (see the proof of Theorem 2). Their result holds for 
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the more general parabolic regions studied in |U] in any dimension. In [2J, 
M. van den Berg used the sharp results of Lifshits and Shi to obtain ana- 
logues of these results for the Dirichlet heat kernel. Finally, similar results 
have been obtained recently by D. DeBlassie and R. Smits [5] for "twisted 
parabolas" in two dimensions. 

In the case of the cone T 7 in R n , the sharp order of integrability of the 
exit position Bt^ is known. For any domain D in W 1 we write B* D for the 
maximum distance of Brownian motion from the origin up to the exit time 
td for D. More precisely, we set 

B* D =sup{|B t |: 0<t <t d }. 

Then, by Burkholder's inequality (Theorem 2.1 in [3]), for any finite, positive 
p there are constants Cz n and C 2 n such that 

(1.3) Cl n E z [nr D + \z\ 2 ] p/2 < E Z {B* T J> < C*„£7, [nr D + \z\ 2 f 2 . 
Thus it follows from (|1.1|) that 

(1.4) E z [B* T J < oo 

if and only if p < p(j, n). We always have \B^ | < Bj, . On the other hand, 
if p > 1 and 

(1.5) E z \B T J<oo 

then by Doob's maximal inequality (1.4) holds, and hence so does (1.1). This 
gives the sharp order of integrability of -E>t 7 in two dimensions. For W 1 , 
n > 3, since T 7 is finite a.s., it is a consequence of Burkholder's Theorem 2.2 
in [3] that E Z [B^ ] p is finite if E z \Bt 7 \ p is finite and this gives the sharp order 
of integrability of Bt^ in all dimensions. The results of Burkholder generated 
considerable interest amongst probabilists and analysts. In particular, M. 
Essen and K. Haliste used harmonic measure techniques to obtain some 
generalizations of Burkholder's results. 

The problem of obtaining the sharp order of integrability of the exit 
position B TV for the parabola V, and of the related random variable B* , 
suggest themselves - a problem that we address in this paper. As we shall 
see below, the order of integrability of B* p and r-p are, as in the case of 
cones, also closely related. We consider, as in jS] and ^0], more general 
regions in R n of the form 

(1.6) V a = {(x, 7)eKx R 11 - 1 : x > 0, \Y\ < Ax a }, 



2 



with < a < 1 and A > 0. (The case a = 1 is the cone for which, as 
described above, we know everything.) We write r a for the exit time from 
V a . 

We begin with a very simple result which shows that the random vari- 
ables B Tp and B* Trp share the same integr ability properties. More precisely, 

Theorem 1. Suppose that a and b are positive constants and that z E V a . 
Then E z [exp[6 \B Ta \ a ]] < oo if and only if E z [exp[6 (B* a ) a }] < oo. 

Naturally, we will need to estimate the distribution functions of B Tv and 
B* for large t in the manner of H1.2|) . The probability that the Brownian 
motion exits a parabola-shaped domain V a outside the ball of radius t, that 
is P z {\B Ta | > t}, is the harmonic measure of that part of the boundary of Pql 
lying outside the ball -6(0, i) of radius t, taken w.r.t. the domain V a . The 
larger quantity P z {B* a > t} is the harmonic measure of the intersection of 
V a with the sphere of radius t taken w.r.t. the intersection of V a with the ball 
B(0, t). In fact, 1 — P z {B* a > t} is the probability of exiting V a without ever 
exiting the ball -B(0, t). These interpretations of the distribution functions of 
B Ta and B* a facilitate the use of some well-known and quite precise estimates 
of harmonic measure. 

Our result for parabola-shaped regions in the plane provides a complete 
solution to the problem and we state it separately. 

Theorem 2. For the exit position B Ta from the parabola- shaped domain V a 
in the plane, 

(1-7) Hm t^\^\P z {\B Ta \ > t}] = - * 

t^oo 2A(1 — a) 

Furthermore, 

(1.8) [expf&lSj 1 -"]] <oo 



if and only if 



b < 



2A(1 - a)' 

Then, E z [exp[6 \B Ta \ a ]] is integrable for each positive b if a < 1 — a, and is 
not integrable for any positive b if a > 1 — a. 

We note that we can determine whether or not exp[6 |-B T( J a ] is integrable in 
all cases, including the critical case a = 1 — a, b = tt/\2A(1 — a)] . 
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For parabola-shaped regions in R", it is proved in [5] that 
(1.9) -B 1 < liminftSr log [P 2 {r Q > t}] 

Q-l 

< lim sup t«+! log [P z {r a > t}] < -B 2 



t^oo 



for two positive constants B\ and B 2 depending on A, on a, and on the 

a-l 

dimension. In [TJ3, the limit of log [P z {r a > t}] is shown to exist and 
its value is determined explicitly. For comparison purposes we observe that 
if n = 2, a = 1/2 and ^4 = 1, the case of the parabola in the plane, the 
results in |10| give 



1 371" 

(1.10) lim t"3 log [P z {t 1/2 > t}] = —z~, 



3vr 2 

' 1 log ./'-It-,.,, > / H - - 

t— »oo 

while it follows from Theorem 2 that 



(1.11) limt-2log P z {\B Tl/2 \ >t} 



-IT. 



We stated Theorem 2 in terms of limits of logs of distributions to draw a 
parallel with the previously cited work on exit times. However, our results 
are sharper than that, as we obtain sharp estimates for the distribution 
itself, (see Proposition 1 below). 

We prove an extension of Theorem 2 to parabola-shaped regions in higher 
dimensions. 

Theorem 3. We let X\ be the smallest eigenvalue for the Dirichlet Laplacian 
in the unit ball o/M™" 1 . For the exit position B Ta from the parabola- shaped 
region V a of 11. 6\) . and for z £ V a , 

a-l , n i . ,11 



(1-12) Hm e-' log [P z {\B Ta \>t}] = — 

t^oo A{1 — a) 

Furthermore, 

E z [exp [ftl^J 1 -"]] 
is finite if b < \/Ai/[^4(l — a)] and is infinite if b > y/Xi/[A(l — a)] . 

Our estimates on the distribution function of B Ta in higher dimensions 
are not sufficiently precise to determine whether or not exp [6 \B Ta | a ] is 
integrable in the critical case a = l-a,b = V^l/ [-4(1 - a)] . This is one 
reason why we state the two dimensional result separately as Theorem 2. 
Moreover, though the method we use to obtain the distribution estimates 
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for B Ta in the planar case are relatively standard in complex analysis, it 
forms the general outline of the method used to obtain lower bounds for the 
distribution in higher dimensions. For this reason also, it seems helpful to 
present the two dimensional case separately. 

A large part of this paper is devoted to adapting the conformal mapping 
techniques that have led to such precise harmonic measure estimates in 
planar domains to the parabola-shaped regions V a in W n . The so-called 
'Carleman method' can be adapted to estimate harmonic measure in these 
domains from above, as we do in Section 4.1. All else being equal, harmonic 
measure is generally largest in the most symmetric case, and so one would 
expect the upper bounds given by the Carleman method to be reasonably 
precise. However, we were not able to use the Carleman method to obtain 
the lower bounds for harmonic measure that we need to determine the exact 
order of exponential integrability of the exit position B Ta . For this, we 
adapt the conformal mapping techniques used to prove Theorem 2. Our 
parabola-shaped domains being symmetric, we can rewrite the distribution 
estimates for B Ta as a distribution estimate in the corresponding planar 
parabola-shaped domain, but at the cost of having to deal with a Bessel-type 
operator rather than the Laplacian. Conformal invariance is lost, but the 
conformal mapping techniques can still be made to work with considerably 
more effort. An outline of our method for obtaining these relatively precise 
estimates from below for the distribution of the exit position in parabola- 
shaped regions can be found at the beginning of Section 4.2. 

Burkholder's inequality 1)1. 3J) . as outlined earlier, allows one to deduce 
the integrability properties of the random variable B?p a for the cone r Q from 
those for the exit time T a for the cone, and vice versa. In the case of 
parabola-shaped domains, we can partially pass from integrability results 
for the exit position to integrability results for the exit time and conversely. 
The proof provides a partial explanation of the connection between the two 
critical exponents, 1 — a in the case of the exit position and (1 — ot)/(l + a) 
in the case of the exit time. 

Theorem 4. We may deduce from 




(1.14) E z exp[6 2 r Q 1 + a ] < oo. 

Conversely, we may deduce from {1.14)/ , with 62 positive, that 11.13)/ holds 
for some positive b\ . 
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2 Proofs of Theorems 1 and 4 

Proof of Theorem 1. Clearly \B Ta \ < B* a and hence if E z [exp[6 (B* a ) a ]\ < 
oo, then E z [exp[6 \B Ta | a ]] < oo, with the same a and b. 

We turn to the proof of the converse. As observed in the Introduction, 
E z t§, < oo, for each finite p. By (|1.3|) we also have E z (B* a ) p < oo for p 
finite. By Doob's maximal inequality, 



E Z {B*J < 



P 



p — 1 



E z \B Ta \ p , forl<p<oo. 



Thus there exists a po such that for all p > po, 
(2-1) E Z (B* T J<AE Z \B T J. 

We choose an integer ko, depending on po and a, such that po < ka for all 
k > k . Then 



E z [exp[b {B*J] 



hi 



uk 00 uk 



ak 



k\ 

k=l 

h) ^k 



k=l 



k=k + l 

00 uk 
k\' 



a k 



k=k +l 

^ l + Y,TM B *J ak + 4^[exp[6|B r J fl ]] 



h) uk 



k=l 



< OO, 

which proves the theorem. 



□ 



Proof of Theorem 4- We argue as in the proof of the classical good- A in- 
equalities, see [3] for example. Suppose Ijl.l4|) holds for some 62 > 0. We 
set (3 = 1 + a. and note that < (3/2 < 1. Then 

P z {B* Ta > t} < P z {B* Ta >t,r a < tf 3 } + P z {r a > t 13 }. 

For the second term Q1.14JI . together with a Chebyshev style argument, gives 



(2.2) 



Pz{T a > t P } < Cexp 



C exp [—62^ 



l-al 
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Take t > 2\z\, so that t — \z\ > t/2. By our assumption on t and by scaling, 



(2.3) P z {B* a > t, r a < t p } < P z { sup \B s -z\ > -,r a <t 

k 0<S<T Q ^ 



t 



< P z < sup \B S — z\ > 

0<s<t' 3 



P < sup \B S \ > I- 

[o<s<t^ 1 

P \ sup \B S \ > -t 
Lo<s<i i 



1-/3/2 



< C exp 



--t 



2-/3 



C exp 



In the last inequality we used the well-known fact (see pj, inequality (2.15)) 
that 

P \ sup \B S \ > \ \ < Cexp [-A 2 /2] . 

It now follows from (|2.2|) and 1)2. 3JI that 

P z {B* Ta >t}<Cexp [-2M 1_a ], 
for all large t, where 2b\ = min{l/8,&2}- Since 

/>oo 

E z [exp^^J 1 -]] = 1 + b / e 6t P,{B* a > t 1 /^} dt, 

Jo 

we obtain 1)1.13)) . 

Conversely, suppose that (|1.13|) holds, so that 

(2.4) P z {B* Ta >t}< Cexp [-hi 1 "**] . 

Let (3 = 1/(1 + a). We assume, as we may, that t is very large and set 

V(t) =V a nB(0,t 13 ) 
and write f for its exit time. Clearly f < r a and 

(2.5) P z {r a > i] < P z {r a >t,r a = f} + P z {r a >t,r a > f }. 

Let us denote the ball centered at and of radius r in R ra_1 by i? n _i(0, r) 
and the exit time of Brownian motion from it by TB n _ 1 (p,r)- We recall that 
for t > 1, 

^{tb^o.i) > t] < Cexp [-Ait] , 



where C is a constant independent of t and Ai is the first Dirichlet eigenvalue 
of £? n _i(0, 1). By our definition of the region V(t), 



(2.6) 



P z {r a >t,r a = f} < P z {f > t} 

< Po{T~B n - 1 (O,t°'0) > *} 

= Po{r Bn _ im > t 1 ^} 



< C exp 
= C exp 



-Ai t 1 - 2 ^ 



-Al tl+a 



where we used scaling for the first equality above. On the other hand, using 

a. 



(2.7) 



^{r a >t,r a >f} < P t {B*>tP} 



< C exp 
= C exp 



. felt /?(l-a) 



The estimates ()2.5[) ~ (j2.7p prove that (|1.14l) may be deduced from (|1.13|) . 
This completes the proof of Theorem 4. □ 

We note that our argument above proves that if E z [exp [62 T a a ]] < 00 for 
some < a < 1 and b 2 > 0, then E z [exp [61 (B* a ) 2a ^ 1+a ^]\ < 00 for some 61 
depending on 62- Conversely, if E z [exp[&i (i?* a ) a ]] < 00 for some a and b±, 
then 

E z 



< 00 



for some 62 depending on 61. Thus by Theorems 2 and 3 we see that 
(2.8) E z [exp[6r Q p ]] < 00 

for some b > if and only 



P — 



1 — a 
T+a' 



as was already proved in |T| for a = 1/2 and in [S] for general a in (0, 1). 



3 Parabola-shaped domains in the plane 

3.1 Harmonic measure estimates 

For the moment we restrict ourselves to two dimensions. The key estimate 
for the distribution function for B Ta , of which the limit (|1.7|) is a direct con- 



S 



sequence, is Proposition 1. As we shall see below, we may suppose without 
loss of generality that z = zq, where zq is the point (1,0). 

Proposition 1. There are constants C\ and C2 depending only on a and 
A such that, as t — > 00, 



C\ exp 



7T 



A— a 



2A(1 - a) 



< P Z0 {\B Ta \>t} 



< C2 exp 



7T 



2A(1 - a 



t L - a + o{t L ~ a ) 



The first step in the proof of Proposition 1 is to show that there is 
a negligible difference between P Zo {\B Ta \ > t} and P Zo {E>l a > t}, where 
for z = (x,Y) € M. x ]R™ _1 we write z 1 for x, the projection onto the 
first coordinate. This will follow from a simple estimate that holds in all 
dimensions. 

Lemma 1. For the exit position B Ta of Brownian motion from V a in W 1 , 
and for sufficiently large t, 

P Z0 {B l Ta >t] < P Z0 {\B Ta \>t} < P ZQ {B l Ta >t-Ah 2a - 1 }. 

Proof. For fixed t positive, we write 5(0, t) for the sphere center and 
radius t and we write x{t) for the common first coordinate of the points of 
intersection of the boundary of V a with 5(0, t). Then 

(3.1) P Z0 {\B Ta \>t} = P Z0 {Bl>x(t)}. 
We also have, for all sufficiently large t, 

(3.2) t - A 2 t 2a ~ x < x(t) < t. 

The upper bound is clear. For the lower bound, we observe that any point 
(x,Y) with x = t-AH 2a - x and \Y\ = A(t - A 2 t 2a ~ l ) a lies inside the sphere 
5(0, t). Indeed, 



+ A 2 {t- A 2 t 2a ~ l ) 2a -t 2 

<4.4a-2 



AH" a ~' z - 2A 2 t 2a + A 2 {t - A 2 t 2a - V ) 2a 

< AH 4a ~ 2 - 2A 2 t 2a + A 2 t 2a 
= A 2 t 2a (A 2 t 2 ^ - 1) 

< 0, 

where we use the assumption that < a < 1 and note that < t—A 2 t 2a ~ x < 
t. □ 
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Since (t - Ah 2 *- 1 ) 1 a = t l ~ a [l + o(l)], Proposition 1 would follow from 
Lemma 1 together with 



(3.3) C\ exp 



TT 



2A(1 - a) 



1-a 



< P ZQ {Bl a >t} 

< C2 exp 



TT 



A—a 



2A{1 - a) 



3.2 Proof of (3.3) by conformal mapping 

We wish to use the Ahlfors-Warschawski estimates on conformal mappings. 
Following the standard notation in this area, we introduce the functions 
tp+(x) = —ip-(x) = Ax a for x > and 6(x) = ip+(x) — ip-(x), so that 
9{x) = 2Ax a for x > 0. Then V a has the form 

T^a = {z = x + iy : x > and <f-(x) < y < ip + (x)}. 

We write f(z) for the conformal mapping of V a onto the standard strip 

S = {w: \lmw\ < tt/2} 

for which f(zo) = and f'( z o) > 0- The mapping / is then symmetric in the 
real axis. We write Et for 9"Po,n{Rez > t} and write F s for SSnjReit; > s}. 
For t > 0, Et corresponds under the mapping / to F s , for some s that we 
denote by s = s(t). 




F s (s = s(t)) 



Then, by conformal invariance of harmonic measure, P Zo {B^. a > t} is the 
harmonic measure at of F s t t \ with respect to S, that is 

(3-4) P Zo {B 1 Ta >t} = co(0,F s(t y,S). 

Harmonic measure in the strip is easy to estimate, for example by mapping 
the strip onto the unit disk where harmonic measure at the origin coincides 
with normalized angular measure on the circle. One may show that 

w(0, F s ; S) = - arg [e 2s - 1 + 2ie s ] , 
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whence, for s large, 

(3.5) -e" s < u(0,F s ;S) < -e~ s . 

7T 7T 

Warschawski's estimates |15| Theorem VII] for the real part of the mapping 
/ involve an error term f°° 6'(x) 2 /9{x) dx for which, in our case, 

yn dx = A i ^dx = 2a 2 ^/ x a_2 dx<oo. 



<9(x) A 2x Q 

From [I3J Theorem VII] we deduce that, for z\ = x\ +iyi and Z2 = X2 + iy2, 

f X2 dx 

Re/(z 2 )-Re/(z 1 )=vr J _ + o(l), 

as sci, X2 — > +00, uniformly with respect to y\ and yi- As a consequence we 
obtain that as t — > 00, 

«•/<«+*>=*/, «w + °w = 2i(r^) tl -° + o( 1 )- 

Since this estimate is uniform in y, it follows that 

(3.6) s ( t ) = _J!_ t i-« + 0(l), asi^oc. 

We will use this estimate on the boundary correspondence under / again in 
Section 4. Together IpP]) . and (|3~6|) yield Thus Proposition 1 is 

proved. 

Note It is possible to shortcut the above explicit calculations by using 
Haliste's estimates for harmonic measure [7J Chapter 1.2], that are them- 
selves based on the Ahlfors-Warschawski approximations. Haliste formu- 
lates her estimates in terms of the harmonic measure of a vertical cross 
cut Ot = [t — iAt a ,t + iAt a ] with respect to the truncated domain V a (t) = 
{(x,y): < x < t, \y\ < Ax a }. We set 

B% = sup{\Bl\: 0<t<r a }. 

Then, 

cu(z ,9 t ;V a (t)) = P Z0 {B^ > t}, 

for which Haliste provides estimates similar to (|3,3|) . A version of Lemma 

1 * 

1 for the maximal functions B T ' a and B* a yields Proposition 1 with \B Ta \ 
replaced by B* . At this point, it suffices follow the argument in the next 
section with B* a instead of \B Ta \ and to recall that \B Ta \ and B* a have the 
same integrability properties (Theorem 1). 
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3.3 Proof of Theorem 2 



As above, we assume for the moment that z = zq = (1,0). First, the limit 
(|1.7|) is a direct consequence of Proposition 1. Next we observe that 

POO 

(3.7) E Z0 [exp[b\B Ta \ a ]] = l + b e bt P ZQ {\B Ta \ > t l ' a }dt. 

Jo 

Proposition 1 yields that 



C\ exp 



7T 



2A{X-a) 



< e bt P Z0 {\B Ta \>t 



< C2 exp 



l/(l-a) 



} 



7T 



2A(1 - a) 



t + o(t) 



We suppose first that b < tt/[2A(1 — a)] and set 

b. 



2A{1 - a) 



Then s is positive and we observe that — et + o(t) < —et/2 for sufficiently 
large t. Hence, for t large, 



e ht P Z0 {\B Ta \>t l l { ~ 1 -°t}<C 2 e- 1 i 



t-t 



It follows from the case a = 1 - a of (J37J) that E Zo [exp[6 \B Ta j 1 " ]] < 00 m 
this case. 

In the case b > n/[2A(l - a)}, 



e bt P Z0 {\B Ta \>t 



l/(l-a) 



} > C\ exp 



7T 



2A{\ - a) 



t 



>Ci. 



This gives E zo [exp[6 |-B Ta | 1_a ]] = 00, for such b. 

The cases 0<a<l — a, 6>0 and a > 1 — a, b > may be handled 
similarly, or one may compare the expected value with that of exp [b \ B Ta 
where, say, b = 7r/[4A(l — a)] and b = 7r/[vl(l — a)], respectively. 

We now remove the assumption that z = (1,0). We first deal with 
the upper bound. We let z = (x, y) and assume that t is very large and 
certainly much larger than 1. By translation of paths it is clear that 
Lo((xi,y),9t;V a (t)) < uj((x2, y),9t]V a {t)) whenever x\ < X2- Hence, we 
may assume that x > 1. Thus for t » x » 1 we have, by symmetry, 



(3.8) 



sup uj(z',0 t ;V Q (t)) = ui(x,6 t ;V a {t)). 

z'&9 x 
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From this and the more general upper bound of Haliste [7j we have that for 
z = (x,y) and t » x » 1, 



oj(z,6{t);V a (t)) 
(3.9) 



< Ci exp 
= C2 exp 



-TT 



du 

, e(u) 

TT 



t 



1-Q 



TT 



„l-a 



2A(1 - a) " 2A(1 - a) 

give the desired upper bound esti- 



These two inequalities, (3.8) and (3.t 
mate on P Z {B\'* > t}. 

For the lower bound, we may assume by the above argument that z = 
(x,y) with < x < 1. As before, we may assume that t is much greater 
than 1. By a standard Whitney chain argument and the Harnack inequality 
we have 



(3.10) uj(z,e t ;V a (t)) > C(z)u(z ,e t ;V a (t)), 

where C{z) is a function of z depending on the distance of z to the boundary 
of Vol- The general lower bound follows from this and the case of z = (1,0) 
which we have already done. □ 



4 Parabola-shaped regions in W 1 

The parabola-shaped regions with which we work have the form 

V a = {(x, 7)eRx R 11 - 1 : x > 0, \Y\ < Ax a }, 

for < a < 1 and A > 0. Our objective is to derive estimates for the 
distribution function of the exit position of Brownian motion from such 
regions or, equivalently, for the harmonic measure of the exterior of the 
ball of center and radius t with respect to such regions. In Section 4.1, 
we derive an upper bound for the distribution function by means of the 
Carleman method and in Section 4.2 we introduce a new conformal mapping 
technique to derive an equally sharp lower bound. 

4.1 Carleman method: Upper bound for harmonic measure 

It will be more convenient to write the domain V a in this section as 

V a = {(x,Y) eRx E"- 1 : x > 0, Y G B n _i(0, Ax a )}, 
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where B n -\(0,r) is the ball in R n_1 centered at and of radius r. For 
convenience of notation we set 6(x) = B n _i(0, Ax a ) and refer to 9(x) as a 
cross cut of P a at x. For t large we set 

V a {t) = {{x, Y) : < x < t, Y e 9(x)}. 

This is the domain V a truncated to the right of t. For such a t and any 
(x,Y) £ V a (t), we denote by ui((x,Y),0(t);V a (t)) the harmonic measure of 
(i, 0) + 6(t) at the point (x,Y) relative to the domain V a (t). We wish to 
estimate w((l, 0), 9{t);V a (t)). 

Proposition 2. There exit two constants C\ and C 2 , that depend on n, X\, 
A and a, such that for t > C\, 



u((i,o),0(ty,r a (t)) <c 2 t< n -^ 2 



exp 



t 



l-a 



A{1 - a) 

Before we begin the proof of Proposition 2, we show how it leads to an 
estimate for the distribution of B Ta . With the notation of Lemma 1 and of 
the note at the end of Section 3.2, 

^(1,0) i\B Ta \ >t} = P m {B l Ta > x(t)} < P m {B% > x(t)}. 

The distribution function for B^'* is the harmonic measure that is estimated 
in Proposition 2, so that 

P m {B%>t} = u((l,0),e(t);V a (t)). 

From the estimate (|3.2|) . namely t — A 2 t 2a ~ l < x(t) < t, it follows, as in the 



proof of Lemma 1, that x(t) 1 a 
2, this leads to 



t a [l + o(l)]. Together with Proposition 



P(i,o){B£ > x(t)} < C 2 x{t)^' exp 



)«(«-l)/2 

< C 2 t a{n -^/ 2 exp 



A(l - a) 



x(t) 



l-a 



A(l-a) 



t i - Q [l + o(l)] 



We may absorb the term C 2 £ Q ( n_1 )/ 2 into the o(l) term in the exponential 
to deduce an estimate for the distribution function for B Ta in the following 
form. 

Proposition 3. Suppose that e is positive. There exists a constant C\ de- 
pending on e, n, \\, A and a such that, for t > C\, 



P(i,o) {\B Ta \ >t}< exp 



A(l-a] 



[l-e]t 



l-a 
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Proof of Proposition 2. Our estimates follow those of Haliste [7] . We take t 
to have some large, fixed value and set 

(4.1) Kx)= uj 2 (x,Y)dY, 0<x<t, 

Je(x) 

where for convenience we write to(x,Y) for u}((x,Y),9(t);V a (t)). Differen- 
tiating h (see Haliste [7] for the justification of this step), we obtain 



(4.2) h'(x)= 2u x (x,Y)uj(x,Y)dY 

Je(x) 

and 

(4.3) h"(x)=2 I oj xx (x,Y)uj(x,Y)dY + 2 / \u x (x, Y)\ 2 dY 

Je(x) Je(x) 

We observe that, since uj(x,Y) is increasing in x for each Y, the derivative 
lo x (x,Y) is non negative and hence h'{x) > 0. Since uj(x,Y) is harmonic, 
we have 

uj xx (x, Y) + Ay uj(x, Y) = 0. 

Thus, 

(4.4) h"{x) = -2 f A Y uj(x,Y)uj(x,Y)dY + 2 [ co x (x,Y) 2 dY. 

Since the harmonic measure vanishes on the lateral boundary of the domain, 
uj(x,Y) = if Y G d9(x), with < x < t. Thus, integrating by parts, we 
obtain 

(4.5) h"(x) = 2 [ \V Y uj(x,Y)\ 2 dY + 2 [ u x (x, Y) 2 dY. 

Je(x) Je(x) 

Writing 5(0, r) for 2? n _i(0, r), we now recall that for all u that are differen- 
tiable on B(0,r) and vanish on dB(0,r), 



(4-6) A B(0 ,r) < 



Jb(Q,t) l Vn ! 2 



Jl 



B(0,r) 



III, 



where Ag(o,r) is the first eigenvalue of B(0,r) for the Laplacian. By scaling, 

As(0,r) = ~2^1j 
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where Ai is the eigenvalue of the unit ball. In our case this gives 

( 4 - 7 ) M*) = Al ' 

From (|4.5[) and (|4.6() we deduce that 



h"(x) > 2\ e(x) / uj{x,Y) 2 dY + 2 / u x {x,Y) 2 dY 
Je(x) Je(x) 



(4.8) = 2\ e[x) h{x) + 2 Lu x (x,Y) 2 dY. 

J6(x) 

On the other hand, by (|4,2|) and Holder's inequality, 

h'{x)<2( f u x (x,Y) 2 dY] 1 ( [ u(x,Y) 2 dY^ 
V Je(x) J V Jo(x) 



or 



h'(x) 2 <4[ [ u x (x,Y) 2 dY)h(x) 
J9(x) J 



or 

h'{x) 2 



d(x) ) 

< 2 f co x (x,Y) 2 dY. 
Je(x) 



2/i(x) Je(x) 
This and ijOjl give 

(4.9) > 2\ e ^h(x) + ^ 

Since 2yfa\fb < a + b, we conclude that 



/i"(x) > 2tJ\ 6 ( x ) h'(x) 
which, by (|4.7j) . is the same as 

h"{x) > 1 
/i'(sc) ~~ A x a ' 

Following Haliste, we consider the function g on the interval (0, t) given by 



g(x) 



J exp ^2 J A e(r) dr J (is 

* /2^AT f dr\ _ 
eXP l^/o ^J* 

exp — — s ds. 
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This function satisfies 



1 J \J X Hr) dr 



g'(x) = exp ^2 

g"(x) = g'{x)2^X^ 
2^\[ 1 



and so 



dx 



(log/i'- logs') > 0. 



From this it follows that the function is non-decreasing on (0, t). Since 
g(0) = h(0) = 0, the generalized mean value theorem gives that for any 
< x < t there is a £ <E (0, x) such that ^ = Hence, 

h(x) h'jx) 
g{x) ~ g'{x) ' 

Since g'(x) > 0, this shows that is nonnegative and hence the func- 

tion ^||y is non-decreasing. Thus, 



(4.10) 



h(x) < h(0 ^ , for < x < e < t. 



Setting //(r) = 2 V / A e(r) , 

<7(x) = exp ^ J /i(r) dr 

(4.11) = exp 

We may estimate H(x) by 
H(x) = 



exp 

J /i(r) dr^ H(x). 



/j,(r) dr I (is 



o 

= exp 
x. 



ex P I - — — [ 



A(l 



A(l-a) 
- ") J Jo 



3.1-a _ g l-a] 



(is 



exp 



2^ 



.4(1 -a) 



1 s l-a 



(is 



Therefore, 

(4.12) (/!./•) < x exp ^ j jiir ) dr ) . < ./■ < / 
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On the other hand, setting 

~A{1 - a) 



i V(i-a) 



x = x Q (\i,A,a) = In 2 + 1 

we find that, for x > xq, 

H(x) = expf-ifx 1 "" 



and K 



A(l 



exp [K 



s 1 -*] ds 



> 



exp f - K x 1 '*] [ — *— r a/ ^- a) exp \K r] dr 
Jo 1 - " 

i i_ r la 

exp I" — K x 1 a ] / exp \K r] dr 

I- a J l 

exp [ - if [ exp [if X 1 ' ] - exp if] 



> 



K(l-a 
1 

if(l - a) 
1 



1 - exp [K (1 -x 1_a )] 



2if (1 - a) 
We have shown that, 



H(x) > 



A 



for x > xq, 



and this, together with ()4.11|) . gives 

(4.13) g{x) > — -= exp [ / /i(r)dr J, for x$ < x < t. 

From (|4TTn|) . dUEl) and (jiTSj) we deduce that, 

M x ) ^ x /i(^) exp ^ — J fi(r) dr^j , for xo < x < £ < t. 

Taking x = x$ and letting £ tend to t, we arrive at 



frOco) < x h{t) exp ^ - J 



ji(r) dr 



C±(\i, A,a) h(t) exp 



A{1 - a) 
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for an appropriate constant C\. By our definition of h, 

h(t) < vol{B n ^ 1 (0,At a ))= ln A n - 1 t^ n - 1 \ 
where 7 n is the volume of the unit ball in Thus, 

h(x ) < C 2 (n,\ 1 ,A,a)t a( - n -V e xp 

for all sufficiently large t. We now choose an r, independent of t, such that 
the ball B n [(xQ, 0), 2r) is contained in V a \ clearly such an r exists. Then, 
by the Harnack inequality, 

w(x ,0) < C 3 (n)u(x,Y) for all (x,Y) G £ n ((x , 0), r) . 

Squaring and then integrating over the ball -B n _i(0,r) leads to 

uj 2 (x ,0) < C 4 (n,r) f uj 2 (x ,Y)dY 

< C 4 (n,r) I u 2 (x ,Y)dY = C 4 (n, r) h(x ) . 
Je(x ) 

From this we finally obtain 

Lo((x ,0),9(t);V a (t)) < C^X^A^)^- 1 ^ 2 exp 

for t > xq. One final application of the Harnack inequality to move from 
(xo,0) to (1,0), and the proposition is proved. □ 

4.2 Conformal mapping method: Lower bound 

In this section we write "P™ for the region in (|1.6|) to emphasize the dimension. 
We observe that, because of the cylindrical symmetry of (this region is 
invariant under rotation about the x-axis) and because of the symmetry of 
the boundary values of the harmonic measure, the value of the harmonic 
measure at (x,Y) in V™ depends only on x and on \Y\. 

We associate with the corresponding domain V a = V 2 in two dimen- 
sions. The technique we develop to obtain lower bounds for the distribution 
function of the exit position in V™ is the following. We replace Laplace's 
equation in T 7 ^ by the corresponding Bessel-type partial differential equa- 
tion in V a ~ in the case n = 2, this reduces to the Laplacian. Mirroring the 



2VAT 
Ail -a) 



tl—a 



A(l 



a 
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arguments in Section 3.2, we map the parabola-shaped domain V a confor- 
mally onto the standard strip S = {w : \Imw\ < vr/2} and determine the 
form the partial differential equation takes in the strip after this change of 
variable. Adapting the Ahlfors-Warschawski conformal mapping estimates 
to our purposes, and in some instances refining them, we show that the 
partial differential equation in the strip is but a small perturbation of the 
Bessel-type partial differential equation that we started with - it is almost 
conformally invariant. The solution of the unperturbed Bessel-type p.d.e. in 
the strip may be easily estimated. We show how to compare the solutions 
of the perturbed p.d.e. with those of the unperturbed p.d.e., so as to obtain 
the estimates we require on the original harmonic measure in the parabola- 
shaped region in W l . As well as keeping track of how the p.d.e. changes as 
we change from one domain to another, we also need to keep track of the 
boundary conditions - but here Warschawski's detailed conformal mapping 
estimates are exactly what we need. 

We break the proof into a number of subsections and lemmas. 

4.2.1 From a parabola-shaped region in IR n to a planar strip 

To begin with we compute how the Laplace operator changes as we drop 
down from n dimensions to 2 dimensions. 

Lemma 2. Suppose that H(x,Y) is a C 2 '-function on that is invariant 
under rotation about the x-axis, so that H depends only on x and on \ Y\. We 
associate with H a function h(z) in the half parabola- shaped planar domain 

'Pa = {z = x + iy: x > and < y < Ax a }, 

defined by 

h(x + iy) = H(x,Y) whenever \Y\=y. 

Then 

(4. 14) AH (x, Y) = Ah(x + iy) + (n - 2) ^ (x + ^ . 



Proof. In this proof, we denote a point in "P™ by (xi, X2, ■ ■ ■ , x n ). With this 
notation, 

H(x 1 , x 2 , ■ ■ ■ , x n ) = h ^xi, \Jx~l + ... + x{^j = h(x, y) 
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with x = x\ and y = \/ x' 2 , + . . . + x\ . Then, 



d 2 



i=l 1 



dx 
dx 2 
dx 2 



i=2 
n 

+£ 

i=2 

dh 
dy 



dh 



1 ■■. 



dy y/xj + ... + x. 



dh {x\ + . . . + x 2 ni 



+ 



d 2 h 



dy {x\ + . . . + x 2 ) 3 / 2 dy 2 x\ + . . . + x 2 

d 2 h 
dy 2 



V4 + 



V4 + 



~\~ x r 



+ 



Ah(z) + (n - 2 



h y {z) 



□ 

From now on we may work in two dimensions and have conformal map- 
ping at our disposal, at the expense of having to deal with the more compli- 
cated Bessel-type differential operator appearing in ()4.14j) . rather than the 
Laplacian. The complication arises because this operator is not conformally 
invariant. 

As in Section 3.2, we denote by w = f(z) the conformal mapping from 
the domain V a onto the standard strip S, for which /(l) = and /'(l) > 0. 
Since / is real on the real axis, the upper half of the parabola-shaped 
domain V a is mapped to the upper half S + of the strip, specifically, S + = 
{w:0<lmw< vr/2}. 

We denote the inverse mapping of f(z) by g(w). We associate a function 
k(w) in S + with a function h(z) in according to 



(4.15) 



k(w) = h(g(w)), for w G 



Then h(z) = k(f(z)), for z G V£ • In the next lemma, we compute how 
the differential operator on the right of (|4.14|) changes under this change of 
variables. 

Lemma 3. Suppose that h(z) is a C 2 -function in the domain and that 
k is defined in the strip S + by \4-15\) . Then, with g(w) = z, 

2(n-2)- 



. . . , ^h v (z) Ak(w) , , Im \k w (uS) I q'(w)] 
(4.16) Ah(z) + - ' ,! 1 ! 



Im [g(w)] 
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Proof. We recall that h(z) = k(f(z)), and use the formulae for change of 
variable in the complex partial derivatives d/dz and d/dz, as explained, for 
example, in Krantz [SJ Section 1.2]. First, 

Ah(z) = (Ak)(f(z))\f(z)\ 2 , zeV a . 

In general, 

d . ( d d 



d(lmz) \dz dz 

whence 

.fdh dh\ 

w = \-d-z-m) 

§- z [Hf(z))]--^[k(f(z))] 

^(m)f(z)-^{f(z))7{z]). 

ow ow 
Since k(w) is real- valued, k w = k w . We therefore obtain, 
h y (z) = -2Im [k w (f(z))f'(z)]. 

Thus, 



Ah(z) + (n - 2) 



h y (z) 

y 

(Ak)(f(z))\f'(z)\ 2 -2(n-2) 



2 ^lm[k w {f(z))f(z)} 



Im [g(w)] 

On substituting to for /(z) and l/g'(w) for f'(z), we obtain ()4.16() . □ 



4.2.2 Asymptotic estimates for the conformal mapping g 

The success of the transformations introduced in the previous subsection 
depends on being able to simplify the expression on the right hand side of 
(|4.16f) , which is essentially an expression for the Laplacian in the domain "P™ 
in R n transformed to the strip S + in the plane. This is achieved using mod- 
ifications of results of Warschawski , which give asymptotic expressions 
for the conformal mapping g and for its derivative. For some of these we 
draw on [I]. We begin with an estimate for the imaginary part of the map- 
ping g. In this context it is helpful to keep in mind that estimates for the 
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real part of the mapping g, that is for the rate of growth of g, are generally 
more difficult. 

We will adopt the more general situation, as considered by Warschawski, 
of a conformal mapping / of a domain of the form 

D = {z : \lmz\ < 4>(Rez)}, 

where 4>(x) is continuous on the real line, onto the strip S. In our case, 
<j)(x) = Ax a , for x positive. Warschawski 's domains are not necessarily 
symmetric, but the symmetric case is sufficiently general for our purposes 
here. Warschawski writes 9(x) = 2<p{x) for the width of the domain D at x. 
We assume that D has boundary inclination at x = oo in that 

0(^2) — ^(^l) = ( x 2 — %i) as x\, x 2 — > 00. 

This is the case if 4>(x) is continuously differentiable and <p'(x) — ► as 
x —* 00, and includes our parabola-shaped domains V a - 

We take the conformal mapping / of D onto S for which / is real and 
/' is positive on the real axis. We denote the inverse mapping of / by 
g, in agreement with earlier notation. We begin with an estimate for the 
derivative of g. 

Lemma A. Warschawski, [15, Theorem X(ii)] For each p with < 
p < 1, the conformal mapping g of the strip S onto D, for which g is real 
and g' is positive on the real axis, satisfies 



(4.17) g'{w) 



-+o(l) 

7T 



e(Reg(w)), 



uniformly as Hew — > 00 in the sub strip S p = {w : |Imui| < p7r/2}. 

In |15l Theorem X(iii)], Warschawski obtains an asymptotic expression 
for Im g(w) as Rew — > 00. Taking advantage of the symmetry of the domain 
D (so that Imiji(w) = when Imw = 0) and adapting Warschawski's proof, 
we prove 

Lemma 4. The conformal mapping g of the strip S onto D, for which g is 
real and g' is positive on the real axis, satisfies 



(4.18) lmg(w) 
as Re w — > 00 with w 6 S. 



- + o(l) 

7T 



6(g(Rew)) Imw, 
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Proof. By symmetry of the mapping g, it is enough to prove ((4.18(1 in the 
case of Imu> positive. First we show that (|4.18|) holds in each sub strip S p , 
(0 < p < 1). By Theorem 11(a) in [To], 

(4.19) lim &igg'(u + iv) = 0, 

u— >oo 

uniformly in v, \v\ < ir/2. Combining ((4. 19(1 with Lemma 1X1 we find that 



Reg'(w) 



- + o(l) 

TT 



9(Reg(w)), 



for w G S p . Then, for w = u + wq G S p , 



lmg(u + ivo) 



lmg(u + ivo) — lmg(u) 

— 1m giu + iv 

ov 

vo Reg'(u + iv\) 
~+o(l) 



[vi G (0, v )] 



V=V\ 



TT 



9(Reg(u + ivi)). 



We assert that, uniformly in v, \v\ < tt/2, 

(4.20) 6(Re g(u + iv)) = 6{g{u)) + o(l), 

(which is why we need not distinguish between 6(g(Rew)) and #(Re g(w)), 
up to o(l)). This then gives the stated expression for \mg(u + ivo), for 
\vq\ < pTr/2. To prove Q4.2UJI we note that, as at the bottom of Page 290 of 
|15j and as a consequence of the assumption that D has boundary inclination 

at x = oo, 

Reg{u + iv) = g(u) + o(l), as u — > oo. 

On using yet again the assumption that D has boundary inclination at 
x = oo, P~271|) follows. 

Finally, we need to show that (|4.18|) holds uniformly in v. Given e small 
and positive, we take p = 1 — e, so that ((4.18(1 holds for w = u + iv, with 
M < (1 — e)vr/2. In particular, 



,7T 



Im# (it + i-(l - e; 



1 



7T 



+ o(l) 



( i- e) ^M + o (%(«))), 



as u — > oo. Thus, using ((4.20(1 . the image of the sub strip 5 P is a region of 
the form 

6(Rez) 



z : |Im z| < (1 — e) 



+ o(6(Re z)) 
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Since g(w) lies outside this region if v > (1 — e)vr/2, we find that 

e(Reg(w)) 



(4.21) 



Img(w) > (1 — e) 



+ o(0(Re s (w))) 



(1 _ £) %M +o(%W) ). 



On the other hand, if v > (1 — e)vr/2 then, simply because lies in Z), 
(4.22) 



g(RegH) %(n)) 
Img(io) < = h o (1). 



Together (j4.21l) and ()4.22|) imply that, for u sufficiently large and v > (1 

e)7r/2, 



g(g(u)) 



7T 



<2c %(«)). 



It follows that (|4.18f) holds uniformly in v, —tt/2 < v < ir/2, with o (1) 
replaced by Ce. Since e may be as small as we please, the proof of (|4.18|) is 
complete. □ 

We wish to remove the restriction in Warschawski's Lemma^that w lies 
in a fixed sub strip of the standard strip S, at least when D is a parabola- 
shaped domain. In this situation, 

0{x) = 2Ax a , for x > 0. 

Proposition 4. We set g to be the conformal mapping of the strip S onto 
V a for which g is real and g' is positive on the real axis. Then the following 
estimate for the derivative of g holds: 



(4.23) 



- + o(l) 

7T 



%(Reu>)) = 2A 



- + o(l) 

7T 



as u, the real part of w, tends to oo, uniformly in Imro. 



Proposition 4 is part of the main result in in which it is shown that the 
function 



Hz) 



exp 



7T 



2(1 - a) 



[1 -(1 + Bi 



is univalent in V a and maps V a onto the interior D of a Dini-smooth curve 
C lying inside the unit circle \z\ = 1, except for z = —1. A smooth Jordan 
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curve C is said to be Dini-smooth if there is an increasing function uj(x), 
that satisfies the Dini-condition 



1 L0(x) 



dx < DO, 



for which the angle @(s) of the tangent to C, considered as a function of 
arclength, satisfies 

\P(s 2 ) - P( Sl )\ <w(s 2 -ai). 

The proof in 0] that h(V a ) is indeed bounded by a Dini-smooth curve, that 
lies inside the unit circle except for z = — 1, is quite involved. Proposition 4 
follows from this in a relatively straightforward manner. We repeat a version 
of the argument here for the reader's convenience. 
We first note that 



7T 



(4.24) h'{z) = 

We consider the mapping 



5 log 



l + z 



h{z), z G TV 



z e A. 



The unit disk A is first mapped onto the strip S, which is mapped by g onto 
the parabola-shaped domain V a and, finally, this is mapped by h onto the 
inside of a Dini-smooth curve C. As explained in [3], we are now in a position 
to apply |13l Theorem 10.2] and may deduce that has a continuous, non- 
zero extension to the closure of the unit disk. In particular, ^'(1) ^ 0, 
and we note that ^(1) = 0. We can derive information on the behaviour 
of the derivative of g from information on the derivative of We write 
w(z) = log[(l + z)/(l — z)], for z £ A. Using the expression (|4.24|) for the 
derivative of h, 



h'{g{w(z)))g'{w{z))w'(z) 

~B[1 + Bg(w(z))} g'(w(z)) 



-■kB- 



1 



g'(w{z)) *(z) 



We let z 



l + z [1 + Bg(w(z))} a 1 
1 from within the unit disk. Then 



1 



-*'(!), 
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which is non-zero. Hence, 



g'{w(z)) 2_ 

[l + Bg{w{z))] a ^ vr J B : 



as z 



1, z G A. 



It is not difficult to see that (Rez)/(1 + -Bz) — > -B 1 as Re z — > oo with 
z £ V a . Substituting z = g(w) in this limit yields 



Reg(w) 



— as Rew ^ oo with w £ S. 



1 + Bg(w) B 

Since the unrestricted limit as z —* 1 within the unit disk corresponds to 
the unrestricted limit as Reu; — > oo within the strip S, it follows from the 
previous two estimates that 



g'(w) 



[Reg(w)Y 
which is (JH13). 



B c 



2A 



7T 



as Rew -> oo, w £ S, 



4.2.3 Asymptotic form of the differential operator 

Armed with the asymptotics for lmg(w) in Lemma |1| and those for g'{w) in 
Proposition^! we are now ready to derive the asymptotics for the differential 
operator 



(4.25) 



Ak(w) 
\g'(w)\< 



2(n 



Im [k w (w) I g' (w)} 
Im [g(w)] 



which acts on C 2 -functions defined in the strip S and arises in Lemma 
We keep in mind that, for a real- valued function k(w) = k(u + iv), 

k v (w) = — 21m (k w (wj). 

First we compute, using Proposition 01 



-21m 



k w (w) 
g'(w) 



-21m 

-2 



k w (w) 



7T 



o{g{u)) 



7T 



(*(«)) 



[1/7T + o(l)]6 (g(u)) 
Im [^H(l + o(l))] 



--k v (w) +Im [o(l)A^, 
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e(g{u)) 

= ^-y^H[i + o(i)]. 

Using Lemma |I] to estimate the imaginary part of g(w), we find that 



Im [k w (w) / g'(w)] 
Im [g(w)] 



k v (w)[l + o(l)] 



2 {g(u)) 

Similarly, and again using Proposition 



{g(u)) [l/ir + o(l)]6(g(u))v 
1X2 *»(% + o(l)]. 



7T 



[l + o(l)]. 



\gi(w)\* e*[g(u)) 
In summary, the differential operator in (|4.25f) becomes 



(4.26) 



7T 



e 2 {g(u)) 



[1 + o(l)] Ak(w) + ( n - 2) [1 + o(l)] 



Summary to date: The work in the foregoing sections has been leading 
up to the following. Suppose that H is harmonic in the domain "P™ in W 1 
and that H(x,Y) is rotationally symmetric about the x-axis. Suppose that 
the function h in the planar domain is constructed from H according to 
h(x + iy) = H(x,Y), with \Y\ = y. Then, by LemmaEJ h satisfies 



Ah(z) + (n - 2) 1 ^- 



0, 



zeP+. 



From h(z) we construct the function k(w) in the strip S + according to 

k(w) = h(g(w)), w G 

where 5 is a symmetric conformal mapping of the strip S onto the parabola- 
shaped domain V a - The partial differential equation satisfied by k is given 
in the next proposition, which follows directly from (|4.26|) and Lemma 01 

Proposition 5. There is a function e(w) in the strip S, with the properties 
that 

(i) e{w) — > as u — > 00, uniformly in v, 
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(ii) whenever the function k(w) arises from a rotationally symmetric har- 
monic function H in V°^, as described above, then k satisfies the partial 
differential equation 

(4.27) Ak(w) + [n + e(w) - 2]^^ = 0, w G S + . 



Remark 1. If the harmonic function H with which we began had lived in 
a cylinder in M. n (of radius tt/2 and with axis along the x-axis), then the 
associated function h would have the standard strip S as its domain of defi- 
nition. The mapping g would be the identity mapping and so k would simply 
satisfy Ak(w) + (n — 2)k v (w)/v = in this case. Proposition^ may be 
thought of as asserting that k behaves asymptotically as if it derived from a 
cylindrical domain. One may also interpret Proposition [3| as asserting that 
while the differential operator Ah(z) + (n — 2)h y (z)/y is not conformally 
invariant in the same way that the Laplacian is, it is asymptotically con- 
formally invariant. The conformal invariance of the Laplacian was used in 
\3.Ji\j in Section 3.2, and Proposition^ is essentially an extension of this to 
higher dimensions. 

4.2.4 Sub solutions and a maximum principle 

We need to determine the boundary conditions satisfied by a function k 
that is constructed, as in Proposition [SJ from the rotationally symmetric 
harmonic function 

H(x,Y) = P iXtY) {\B Ta \>t}, 

in the region "P™ in M. n . Thus H is the harmonic measure of the exterior 
of the ball of radius t w.r.t. V™. The gradient of H w.r.t. Y vanishes when 
x = because of the rotational symmetry. This translates into the boundary 
condition h y (x,0) = for the associated function h in the planar domain 
V£ and, in turn, into the boundary condition k v (u, 0) = for the function 
k{u,v) in the strip S + : 

(4.28) k v (u) = 0, -oo<u<oo. 

The boundary values of H lead to the condition on the boundary of that 
h{x,Ax a ) = 1 if \(x,Ax a )\ > t and h(x, Ax a ) = if \{x, Ax a )\ < t. Under 
the conformal mapping / of the domain V a onto the strip S, this becomes 
the following boundary condition for k: 

(4.29) k(u + in/2) = 0, -oo < u < s; 
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(4.30) 



k(u + iir/2) = 1, s < u < oo. 



Here s depends on t, as specified in ()3.6|) . We note that the point (1,0) in 
corresponds to the point on the boundary of the strip S + . 
Let us therefore suppose that k is a solution of the p.d.e. 

(4.31) Afc(w) + [n + eW-2]^ = 0, w G S + , 

where e(w) — > as u — > oo, uniformly in v, with the boundary conditions 
(|4^|) and (|43U)) . We will show in the next section that A:(0) decays 
at a slower exponential rate as s — > oo than solutions of 

(4.32) Ak(w) + [n + 5-2]^M = , w£S + , 

v 

when 5 is positive, the boundary conditions being the same as those satisfied 
by k. 

This comparison between the solutions of (|4.31j) and (|4.32|) has a heuristic 
interpretation that may be helpful to keep in mind. In the limiting case 
5 = 0, the differential equation (|4.32|) becomes Ak(w) + (n — 2)k v (w)/v = 0, 
the solutions of which, with the above boundary conditions, may be thought 
of as deriving from harmonic measure in a cylinder of radius tt/2 in M n . The 
solutions of (|4.32|) may then be thought of as corresponding to harmonic 
measure in such a cylinder in a slightly higher 'dimension' when 5 is positive. 
Thus, our results will show that the distribution function of the exit position 
of Brownian motion from a parabola-shaped region in M. n decays like the 
distribution function of the exit position from a cylinder in R n , but with a 
time change that is given explicitly by ()3.6|) . 

It is natural to consider solutions of (|4.32j) in the half strip 

= S + n {u < s} = {w = u + iv : u < s and < v < vr/2}. 

In fact, by symmetry, a solution of the p.d.e. ()4.32|) in S + that satisfies the 
boundary conditions (J4.28|) . (|4.29|) and ()4.30|) will take the constant value 
1/2 on the vertical cross cut u = s. This boundary condition can be satisfied 
by using separation of variables to solve (|4,32|) in the half strip and then 
taking a series expansion. The rate of decay of the solution at as s becomes 
large is then determined by the first term in this Bessel series. This is the 
term that is therefore of interest to us. For each m, we write J m (v) for the 
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Bessel function of order m, we write j m for its smallest positive zero and we 
set 

J m (v) = v~ m J m (v). 

The first term in the Bessel series for a solution of (j4.32|) in is (a constant 
times) 

(4.33) <j> s {w) = e 2jm{u - s)/7T J m v\ , where m = ~(n + 5 - 3). 

Since J m satisfies the differential equation (see Section 17.22], for ex- 
ample) , 



Z(v) + [2m + + J m (v) = 0, 

v 

w) satisfies the p.d.e. (|4.32j) in Sf, as well as the boundary conditions 
Q4.28|) and (|4.29|) . On the vertical side of the half strip Sf, its values are 
simply J m [2j m v /n) . One needs to take the entire series to have a solution 
which equals 1/2 on the vertical cross cut u = s. 
We write L for the operator 

(4.34) L[/]=A/+[n + e(u;)-2]^. 

v 

Then k(w) is a solution of L[k] =0 in the half strip with the boundary 
conditions (|4.28|) and (|4.29|) . In order to compare k to solutions of (|4.32|) . 
we construct sub solutions for L in and then use a maximum principle. 
Of course, all our estimates need to be uniform in s. We show how to obtain 
the sub solutions that we need in the next lemma. 

Lemma 5. We suppose that, for a fixed positive 5, the number us is chosen 
so large that 2|e(u>)| < 5 foru > u$. We suppose that a function k$ is defined 
in the rectangle R s = n {u > us} and satisfies the partial differential 
equation \4-^°^j) there. Suppose further that dks/dv is negative in R s . Then, 
L[ks] > in R s . In particular, L[<ps] > 0. 

Proof. With kg as in the statement of the lemma, 

L[k s ] = Aks+[n + e(w)-2}-^- 

= Aks+[n + 5-2] 1 -^+[e( W )-8] 1 -^ 

V OV V ov 

r r \ ^ ldk& 
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Since 2|e(tt))| < 5 for u > us, it follows that e(w) — 5 has the same sign as —5 
in the rectangle R s . Since dkg/dv is negative in R s , we deduce that L[kg] 
has the same sign as 5 in R s . 

The statement about (f>s now follows from the facts that 4>$ satisfies (|4.32j) 
and that J m is decreasing on the interval (0, j m ). 

□ 

The other ingredient we need is an appropriate form of the maximum 
principle. While the version presented here is most probably not new, we 
have been unable to find it in the literature. Consequently, we outline the 
proof for completeness. 

Lemma 6. If f is a sub solution of L which is C 2 in the closure of the 
rectangle R s (the second derivatives are continuous up to the boundary) and 
which is non positive on the top, left and right parts of the boundary, then 
f(u,v) < for any (u,v) G R s . 

Proof. Let Z% = (Xt,Yt) be the diffusion associated with the operator L. 
Then Yj > almost surely for all t. This is true in the case e(w) = 0, since 
Zt is then a Bessel process and as such it never hits zero (see |1 1| Chapter 
XI]). If e is not zero, we still assume that —5 < e(x,y) < 5 in the rectangle. 
It follows by a stochastic comparison theorem argument (as in the classical 
Ikeda-Watanabe theorem, |12j ) that Yt > almost surely for all t > 0. Now, 
let t be the first time that Zt hits the boundary of the rectangle with the 
diffusion starting at zq = (xo,yo) G R s . This time is finite almost surely. 
Of course, by the above, Z T belongs only to the three sides of the rectangle 
with probability 1. Applying Ito's Lemma, 

(4.35) f(Z tm - mr ) - f(z ) = M t + A t 

where M t is a martingale and A t = J * L[f](Z s )ds. Since / is a sub solution 
of L, we have L[f] > 0. Taking expectations of both sides of (|4.35j) . we 
conclude that 

f(z ) < E Z0 (f(Z tminT )) . 
We now let t — > oo. Since / is bounded in the closure of the rectangle, 

f(z )<E Z0 (f(Z T )), 

which proves f(zo) < 0. □ 
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4.2.5 Rate of exponential decay 

We now have all the ingredients necessary to prove the following estimate 
of k. 

Proposition 6. Suppose that the function e(w), w £ S, satisfies e(w) —* 
as u — > oo ; uniformly in v. Suppose that k(w) is the solution of j4-27\ ) in 
S + with the boundary conditions \4-28{l , \4-29^ and \4-30[) , so that k derives 
from harmonic measure in as in Section 4- £-3. Then, given e positive, 



k(0) > exp 



2jr. 



7T 



+ e 



for all sufficiently large s, where m = (n — 3)/2. 

Proof. The function k is bounded by 1 on the vertical side u = s of R s 
(since it is but a certain harmonic measure in the parabola-shaped region 
in disguise). More precisely, k(s,v) — > 1/2 as s — * oo, uniformly for 
v in (0, 7r/2). In fact, the harmonic measure of dV^ H {x > t} w.r.t. the 
parabola-shaped region approaches 1/2 uniformly on the cross section 
{(t,Y) : Y G M n ~\ \Y\ < At a }. Hence, for all sufficiently large s, 



(4.36) 



- < k(s, v) < 1, for < v < it/2. 



Given e positive, we set m\ = |(n + S — 3) and choose <5 positive, but so 
small that j mi < j m + 7re/4. This is possible since the first positive zero of 
the Bessel function depends continuously on the order of the Bessel function 
and increases with this order. We suppose that is as in Lemma 03 and 
that s > ug. Direct comparison of k with the function <f>s of (|4.33j) doesn't 
quite work, as 4>s is positive on the side u = us of R s while we do not know 
k there. We consider the positive function 



ks(w) 



3 2j mi («-s)/7T _ 2j rni [(u s -s)-(u-u s )]/-ir 



2 J mi 



TT 



w G R s 



in which the second exponential term compensates for the positive values 
of 4>s on the side u = u$. This function is a solution of (|4.32j) in R s and, 
moreover, dks/dv is negative in R s since J mi is decreasing on (0, j mi ). It 
follows from LemmaElthat L[kg] > 0. The function ks satisfies zero Dirichlet 
boundary conditions on the sides u = ug and v = tt/2 of R s , and zero 
Neumann condition on the side v = 0. On the right side of the rectangle its 
boundary values satisfy 



k s (s + iv) 



1 



4jmi(us-s) 



Jmi (2j mi v/ir) < J mi (2j mi v/7r) < J mi (0), 
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for all sufficiently large s. Together with (|4.36|) . we see that we can choose 
a fixed small, but positive, b\ such that 



hk s (s + iv) < 6iJ mi (0) < 1/4 < k(s + iv), 

for all sufficiently large s. Thus, b\ks — k < on the three sides u = us, 
u = s, v = tt/2 of R s , while b\k$ — k satisfies a zero Neumann condition on 
the side v = 0. Since L \b\ks — k] = biL[kg] > in R s , we conclude from 
the Maximum Principle, Lemma |SJ that b\kg — k < in R s . This leads to a 
lower bound for k(u$ + 1) since 

k(us + 1) > hks(us + 1) 

= 2he 2jm i us ^ sinh(2i mi /vr)e~ 2 ^i s /V mi (0) 
> 6 2 e~ 2j ™i s/7r . 

By the Harnack inequality, k(0) > b^k{us + 1), for a constant that does not 
depend on s. Setting 64 = 6362, 



k{0) > b 3 k(u s + 1) > 6 4 exp 



2jm ^ 

S S 

7T 2 



As 64 does not depend on s, we have 64 > e eS//2 for all sufficiently large s. 
The proof of Proposition 6 is complete. 

□ 



4.2.6 Lower bound for harmonic measure 

The point in the strip S corresponds to the point (1,0) in the region 
under the transformations in Subsection 4.2.1. Thus, 

P m {\B Ta \>t} = k(0), 

where the function k satisfies the partial differential equation 1)4. 27|) of 
Proposition 5 and the boundary conditions 1)4.28)) . ()4.29j) and ()4.30)) with 
s = s(t) as given by ()3.6)) . Thus Proposition 6 leads directly to the following 
lower bound for harmonic measure in the parabola-shaped region "P™. 

Proposition 7. Suppose that e is positive. There exists a constant C2 de- 
pending on e, n, A and a such that, for t > C2, 

P(l,0){\Br a \ >t}> exp 



A(l 



a 



[1 + 6] & 
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4.3 Concluding remarks 

The distributional inequalities in Propositions 3 and 7 lead immediately to 
the limit ()1.12j) . and to Theorem 3 by following the steps in the proof of 
Theorem 2 in Section 3.3. 

It is natural to hope that the machinery constructed in Section 4.2 would 
lead to an upper bound for harmonic measure, and not only to a lower 
bound, thus rendering the use of the Carleman method and Section 4.1 
unnecessary. In fact, there is no mention in Section 4.2 of bounds of any kind 
until Section 4.2.5. However, we have been unable to prove a counterpart 
for Proposition 6 involving an upper bound for k(0). 

Acknowledgment. We are very grateful to an anonymous referee for the 
careful reading of this paper and the many useful comments. 
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